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We investigate the impact of the tau and bottom Yukawa couplings on the transport dynamics 
for electroweak baryogenesis in supersymmetric extensions of the Standard Model. Although it 
has generally been assumed in the literature that all Yukawa interactions except those involving 
the top quark are negligible, we find that the tau and bottom Yukawa interaction rates are too fast 
to be neglected. We identify an illustrative "lepton-mediated electroweak baryogenesis" scenario 
in which the baryon asymmetry is induced mainly through the presence of a left-handed leptonic 
charge. We derive analytic formulae for the computation of the baryon asymmetry that, in light of 
these effects, are qualitatively different from those in the established literature. In this scenario, for 
fixed CP-violating phases, the baryon asymmetry has opposite sign compared to that calculated using 
established formulae. 



I. INTRODUCTION 

Electroweak baryogenesis (EWB) is an attractive and testable explanation for the origin of the baryon 
asymmetry of the universe (BAU). Characterized by the baryon-number-to-entropy -ratio ns/s, the BAU has 
been measured through studies of Big Bang Nucleosynthesis (BBN) and the cosmic microwave background 
(CMB) to be in the following range 

{(6.7 - 9.2) X 10"" BBN 
(1) 
(8.36 - 9.32) X 10-" CMB 

at 95% C.L. ||ll|2|]. Assuming that the universe was matter-antimatter symmetric at some initial time (e.g. 
at the end of inflation), the creation of the BAU requires three conditions (the Sakharov conditions yl]): 
(1) violation of baryon number, (2) violation of C and CP, and (3) either a departure from equilibrium or a 
violation of CPT. 
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In EWB, these conditions are realized in the following way. First, a departure from equiUbrium ispro- 
vided by a strongly first order electroweak phase transition (EWPT) at temperature T ~ 100 GeV 
During the EWPT, bubbles of broken electroweak symmetry nucleate and expand in a background of unbro- 
ken symmetry, filling the universe to complete the phase transition. Second, CP-violation may arise from 
complex phases. These phases induce CP-violating interactions at the walls of the expanding bubbles, where 
the Higgs vacuum expectation value (vev) is time-dependent, leading the production of a CP-asymmetric 
charge density. This is the so-called CP-violating source. This CP-asymmetry, created for one species, dif- 
fuses ahead of the advancing bubble and is converted into other species through inelastic interactions in the 
plasma; in particular, some fraction is converted into left-handed fermion charge density, denoted ul. Third, 
baryon number is violated by non-perturbative S'C/(2)7;, processes, which are unsuppressed outside the bub- 



bles, in regions of unbroken electroweak symmetry |^ 
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Following the common usage, we will 



refer to these as sphaleron processes. The presence of non-zero n l biases the sphaleron processes, resulting 
in the production of a baryon asymmetry [11]. Electroweak sphalerons become quenched once electroweak 
symmetry is broken, as long as the EWPT is strongly first order; therefore, the baryon asymmetry becomes 
frozen in once it is captured inside the expanding bubbles. 

In this work, we consider the charge transport dynamics during the EWPT: that is, how charge densities, 
induced by CP- violating sources, diffuse, interact, and get converted into ni, ultimately inducing ri p/s . 



Although in the Standard Model (SM) this dynamics are insufficient to produce the observed BAU 111211 . 
supersymmetric extensions of the SM can readily include all the ingredients to make it successful^. The 
most commonly accepted supersymmetic scenario is the following: the expanding bubble wall leads to a CP- 
violating source for charge density in the Higgs sector, which is then converted into third generation quarks 
through top Yukawa interactions, which in turn is converted into quark charge density of all generations 
through strong sphaleron processes 1 14]. The rate for baryon number production is proportional to rii; in 
this picture, receives contributions from left-handed quarks of all three generations. 



However, as reported in a previous publication jlST . we have observed that bottom and tau Yukawa 



interactions, shown in Fig.[T] cannot in general be neglected from the computation of in supersymmetric 
EWB scenarios. While bottom and tau Yukawa couplings are small in the SM, in supersymmetric extensions 
they can be larger when the ratio of the vevs of the Higgs doublets, Vu/vd = tan /?, is greater than unity. 
The inclusion of bottom and tau Yukawa interactions can change the EWB picture dramatically. 

• Quark charge density may be supressed: For tan /3 > 5, bottom Yukawa interactions become non- 



' Another reason that EWB is not viable in the SM is that there is no EWPT; for a Higgs mass mn > 114 GeV, electroweak 
symmetry is broken through a continuous crossover |4^. 
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FIG. 1: Examples of bottom and tau Yukawa interactions from (a) absorption/decay, and (b) scattering processes 
involving an addition gauge boson, showing how Higgs density is converted into left-handed quark and lepton density. 
Third generation Yukawa couplings also give rise to F-term-induced scattering processes (c). 

negligible, leading to two important effects: (1) strong sphaleron processes no longer induce charge 
densities of first and second generation quarks, and (2) the third generation left-handed quark charge 
density vanishes when the masses of right-handed bottom and top squarks are equal, or when their 
masses are large compared to the temperature T. 

• Lepton charge density generated: For tan (5 > 20, tau Yukawa interactions also are non-negligible, 
leading to the conversion of Higgs charge density into third generation lepton charge density^. 

These novel effects, which come into play for moderate tan [3, can lead to qualitatively different situations 
from those previously considered for supersymmetric EWB. In the present paper, we focus on a new sce- 
nario, where ul can be purely leptonic. We call this "lepton-mediated electroweak baryogenesis." This 
scenario occurs when tan f3 > 20, and the right-handed bottom and top squark masses are, either, approxi- 
mately equal ("t.^^ ~ ^tR^' large compared to T {m-^^ , m-jr^ > 500 GeV). (In other regions of parameter 
space where the quark contribution to is not quenched, the lepton contribution may still provide an ad- 
ditional enhancement or suppression to the total n/,.) 

In the lepton-mediated EWB scenario, the value of ur/s has opposite sign compared to the value of 
ur/s computed when neglecting the bottom and tau Yukawa rates. The ingredients of this scenario will 
be tested in the near future at the Large Hadron Collider and by precision electric dipole moment (EDM) 



searches 



Clearly, to the extent that these experiments can determine the supersymmetric spectrum, and 
the signs and magnitudes of relevant CP-violating phases, inclusion of these Yukawa rates may be essential 
for testing the consistency of supersymmetric EWB with observation. 

In Sec. ini we present the system of Boltzmann equations, generalized from previous work 



Measurements of the muon anomalous magnetic moment favors large tan /3; see, e.g., Ref. |45 
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include bottom and tau Yukawa interactions. In Sec. |TIT1 we provide an analytic estimate of the baryon 
asymmetry in detail. We solve the Boltzmann equations analytically in the limit that tan/? » 1, such that 
bottom and tau Yukawa interactions are in chemical equilibrium. A new qualitative feature of our analysis 
is our treatment of lepton diffusion; we argue analytically how the left-handed lepton charge density (and 
therefore n^/s) is enhanced by virtue of right-handed leptons diffusing more efficiently in the plasma. 

In Sec. |lVl we verify these conclusions numerically. First, we calculate the bottom and tau Yukawa 
interaction rates, showing in what regimes they ai^e sufficiently fast to induce chemical equilibrium. Next, 
after defining the parameters of our lepton-mediated EWB scenario, we solve the system of Boltzmann 
equations numerically. We illustrate all of the aforementioned new effects and verify the agreement between 
our numerical and analytical solutions. In Sec. |Vl we summarize our results. The Appendix summarizes 
some additional numerical inputs used for this work. 



n. BOLTZMANN EQUATIONS 
A. Preliminaries 



The transport dynamics leading to CP-asymmetric charge densities during the EWPT are governed by a 
system of Boltzmann equations. These Boltzman equations have been derived using the closed-time-path 
formulation of non-equilibrium quantum field theory 1 19], leading to a system of equations of the form 

(2) 



ii = — ^ ^^x[lJ-i + IJ-j + ••• - /Xfc - - ••• ) + s\ 



X 



where is the charge current density of the species i. The density j^^^, induced by CP-violating source S^^^, 
is coupled to other species via coefficients Tx that describe the rate for a process i + j + ... k + I + ... 
to occur. (We have explicitly factored out of Tx, for reasons that will become clear below.) The 
chemical potentials are denoted by jii. Chemical equilibrium, occuring when 



Mi + /^J + ••• - l^k- IJ'l 



0, 



(3) 



is maintained when when the interaction rate Tx is sufficiently large 



Following previous work 11141 



ITO, we simplify Eq. (|2]l in three ways. First, we assume a planar 



bubble wall profile, so that all charge densities are functions only of z, the displacement from the moving 



21 



22L 



2311 . which allows us to replace 



bubble wall in its rest frame. Second, we apply Fick's law 
j i —DiVrii on the LHS of Eq. Q, with charge density = j-^. The diffusion constant Di is the mean 
free path of particle i in the plasma. Third, the chemical potentials appearing in Eq. Q are related to their 
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corresponding charge densities by 



ni = —ki^Xi + 0i^^) , (4) 



where we have performed an expansion assuming /Xj/T <^ 1. In the above, the statistical weight ki is 
defined by 
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oo 



in which gi counts the number of internal degrees of freedom, the + (— ) sign is taken for fermions (bosons), 
and the mass of the ith particle is taken to be the effective mass at temperature T. In our analysis to 
follow, these A;-f actors are ubiquitous; they essentially count the degrees of freedom of a species in the 
plasma, weighted by a Boltzmann suppression. 

Through these three simplifications, the Boltzmann equations become a system of coupled, second order, 
ordinary differential equations for the set of charge densities ni{z). Ultimately, it is the total left-handed 
fermionic charge density 

3 

nL = Yl (n^l + + + (6) 

i=l 

that biases weak sphaleron transitions, thereby determining nsj s. 

While in principle there is an interaction coefficient Fx for every interaction in the MSSM Lagrangian, 
we can determine which ones need to be taken into account for the computation of n/, by considering the 
relevant time scales. After a time t, charge densities created at the bubble wall will have diffused on average 
a distance ddiff = \/ T)t (with the effective diffusion constant T) to be defined below). At the same time, 
the moving bubble wall advances a distance dwaii = t. The diffusion time scale, defined by ddiff = c^waih 
gives the time that it takes for charge, having been created at the bubble wall and having diffused into the 
unbroken phase, to be recaptured by the advancing bubble wall and be quenched through CP-conserving 
scattering within the phase of broken electroweak symmetry. This time scale is 

Tdiff = D/vl . (7) 

Numerically, we have Tdiff ~ lO'^/T (shown in Sec. HVl). To this, we compare tx = T^^, the interaction 
time scale. If tx ^ T^m^ then the process f+j + ...<->A; + £ + ...is slow and Tx may be neglected from 
the Boltzmann equations. Physically speaking, charge density is recaptured by the advancing bubble wall 
before conversion processes can occur. On the other hand, if tx ^ T^m, then these interactions are rapidly 
occuring as the charge density is diffusing ahead of the advancing wall, leading to chemical equilibrium ([3]). 
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Expressed in terms of charge densities, the chemical equilibrium condition is 

In this case, the interaction Tx must be included in the Boltzmann equations. 

A similar argument tells us how we expect deviations from Eq. ([Hi to arise. Suppose that species i is 



produced from the expanding bubble wall at z = 0. On distance scales \z\ < \JDtx, close to the bubble 
wall, Eq. ([8]l will break down: particles i have not had enough time to interact via V x- 

B. Setting up the Boltzmann equations 

We now derive the Boltzmann equations within the context of the Minimal Supersymmetric Standard 
Model (MSSM). In principle, the complete system of Boltzmann equations encompasses one equation for 
each species of particle. However, the assumption that certain interactions Vx are in chemical equilibrium 
(such that Tdiff » Tx) implies relations among the relevant chemical potentials (and therefore among 
their corresponding charge densities), allowing one to reduce the system. First, we assume that weak 
interactions (neglecting flavor mixing) are in chemical equilibrium, so that particles in the same isodoublet 
have equal chemical potential. Second, we assume that gaugino interactions (involving SM particles and 
their superpartners) are also in chemical equilibrium, so that a particle and its superpartner have equal 



chemical potential ll29ll . 



Under these assumptions, the complete set of charge densities relevant for the computation of n^/s is 

Di=n^^+n^^, H = n^+ +nHO-n^- -nfjo+ng^+rifjo , (9) 

Ri = n i + , L{ = n i -\- n i + n~i -\- n~i , 

^R '^L '^L 

where i G {1,2,3} labels the generations. Furthermore, we define the following additional notation: 
Q = Q3,T = U3,B = D-i, L = Li, and R = R^. 

The system of Boltzmann equations contains, in principle, a coefficient Fx for every interaction in 
the MSSM. However, interactions that satisfy tx ^ t^is may be neglected. In particular, we neglect 
interactions induced by first and second generation quark and lepton Yukawa couplings. The weak sphaleron 
rate F^s may also be neglected, since t^s ~ 10^ /T S> Tdiff 1361]. Therefore, baryon and lepton number are 
conserved in the collision terms of the Boltzmann equations. 

Not all of the densities in Eq. Q are independent. Neglecting electroweak sphalerons from the Boltz- 
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mann equations, baryon and lepton number are individually conserved: 

3 /TV") 3 



/oo poo 
dz {Qi + Ui + Di)= dzY^ {Li + Ri) = 
-oo .^^^ J-oo -^^ 



(10) 



Because the left- and right-handed (s)lepton have different gauge quantum numbers, they have different 
diffusion constants in the plasma. Even though lepton number is globally conserved, regions of net lepton 
number can develop since R diffuses more easily than L since right-handed (s)leptons do not undergo SU(2) 
gauge interactions. For quarks and squarks, this does not occur since the left- and right-handed (s)quark 
diffusion constants, dominated by strong interactions, are approximately equal lisj]. Therefore, baryon 
number is locally conserved: 

3 

Y,iQi + U^ + D^)=0. (11) 

1=1 

Other simplifications arise since we neglect first and second generation Yukawa couplings. There is 
no production of first and second generation lepton charge, so Li^2 = Ri,2 = 0. Next, first and sec- 
ond generation quark charge can only be produced through strong sphaleron processes, e.g., titu 
^LbR 2 ^i?' changing the number of left- and right-handed quarks by one unit per flavor. 

Since first and second generation quarks are produced in equal numbers, we have 

Qi = Q2 = -2Ui = -2U2 = -2Di = -2D2. (12) 

Together, Eqs. (I11I12I ) imply that 

B = -{T + Q). (13) 

Therefore, we may consider a reduced set of Boltzmann equations involving only the densities Q, T, Qi, 
L , R, H; the remaining densities are then determined by Eqs. (I12I13I ). The equations are of the form of 
Eq. dU, where we use the relation given in Eq. ^ to express the chemical potentials in terms of charge 
densities. For the quarks and squarks, we obtain 
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r 



^h/ V^Q kT/ * 
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kg kr k. 



Qi - Dq Q'{ 



-2Tss (2-^ - — + 



T + Q 1 



i=l 
2 



kq kT ks 
and for Higgs bosons and Higgsinos we have 



i=l 



Qi 
Qi 



Vu, H' - Dh H" 



yt 



— - — + — ^ - r — + s^-'' 

kg kx ku) ku ^ 



^ kn ks kn 



'-^^ '^W kn kn) 



and lasdy for leptons and sleptons we have 

v^L'-DlL" -- 
Vy.B' -DrR" = T 



'^^ \kL kji kuj 



- yr 



f L R 

\kL 

L R H 



kr. kn kn J ^ V kb 



L R 

kh kR 
L R 



5r 



kR 



The relevant interaction coefficients in Eqs. (|14H16I) are as follows: 



(14a) 



(14b) 



(14c) 



(15) 



(16a) 
(16b) 



The coefficients Tyi, where i £ {t, b, r}, denote the interaction rates arising from third generation 
Yukawa couphngs yi. (The top Yukawa interaction rate has been denoted Ty in previous work.) 



The strong sphaleron rate is Tgg — 16 k' ctg T, where oig is the strong coupling and k' ~ ^(1) 11341] . 



• The coefficients F/j and Fmi, where i G {t, b, r}, denote the CP-conserving scattering rates of parti- 
cles with the background Higgs field within the bubble 1 16(1 . 

We also allow for new CP-violating sources S~^_^, although in the present work we do not evaluate their 
magnitudes. In the MSSM, the most viable CP- violating source is arising from CP- violating Higgsino- 



Bino mixing within the expanding bubble wall [371]; in our work, we take this as the sole source of CP 



violation. The constant Vw — 0.05 is the velocity of the expanding bubble wall 1331]. The /c-factors, e.g. 



kR — kru + k^j^ , kg — ktj^ + kbj^ + /c^^ + k^^ , 



(17) 
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follow the same notation as in Eqs. (O. 

After solving the system of Boltzmann equations (|14H16I) for each density, the left-handed fermion 
charge density is 



riL 



where kg = kf^ + fcf,^ , k£ = ky^^ + /c^-i, > ^^c. The three terms in Eq. dTS] ) correspond to the contributions to 
ni from third generation quarks, first/second generation quarks, and third generation leptons, respectively. 
If the masses of all left-handed squarks and sleptons are much above the temperature of the phase transition, 
only fermions contribute to the left-handed density and we have 

nL^Q + 2Qi+L . (19) 

Finally, we show how our Boltzmann equations reproduce those given in previous work in the limit 
Vh: Vt — > 0. In this limit, we can neglect the rates Tyh,T and rmf),r> and CP-violating sources 5~!y. First, 
since there is no source for lepton charge, we have L = R = 0. Second, the only source for B density is 
strong sphaleron processes; therefore, we have 

-2B = Qi (20) 

in analogy with Eq. Thus, Eqs. (1 13 1201 ) imply that Qi = 2{Q + T). Therefore, by Eq. ([H), we have 
the often-used relation rii = 5Q + 4T; this relation is no longer valid for Tyi,,Tyr < Tdiff. In addition. 



the Boltzmann equations of Refs 

TybjTyr < Tdiflf. 



171] follow from Eqs. (ll4a[ bJT5]): they too are no longer valid for 



III. ANALYTIC RESULTS 

In this section, we estimate the solution to the Boltzmann equations, Eqs. (I14H16I ). of which the endpoint 
is an expression for ul, the left-handed fermion charge density that biases weak sphalerons. We assume that 
top, bottom, and tau Yukawa interactions, in addition to strong sphaleron and gaugino interactions, are all 
in chemical equilibrium. These assumptions lead to a series of conditions relating the chemical potentials, 
and therefore the number densities, of various species. By exploiting these relations, we will express all 
quark and lepton densities Q, T, Qi, L, R in terms of the Higgs density H; then, we will simplify the fuU 
system of Boltzmann equations to a single equation for H, which is analytically solvable lil4ll . 
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A. Lepton charge densities 

When tau Yukawa interactions are in chemical equilibrium condition, the relation 

ki ku kR ^ ' ^^^^ 
is satisfied. The sum of the Boltzmann equations for L and R ( [T6l ) is 

Vu,{R + L)' -{DrR" + DlL")=^. (22) 

Since the left- and right-handed lepton diffusion constants are not equal, there is no simple relation that 
would allow us to relate R to L. However, in the static limit (where Vw —>■ 0), Eq. (l22l ) implies that 

DlL = -DrR. (23) 

(We have assumed the boundary conditions L(oo) = L'(oo) = R{oo) = R'{oo) = 0.) Therefore, we have 

L{z) = H{z) + AL(z) = ^ ^l^^ H{z) + AL{z) (24a) 

KH Dl ki + Dr kR 

R{z) = KR H{z) + AR{z) = - ^ /^^^^ H{z) + AR{z) , (24b) 

kn Dl ki + Dr kR 

where AL and AR are the corrections to these relations, derived below. 

Let us now describe the physics of Eqs. (l24l ) through two limiting cases. Case (i): set Dr = Dl- In this 
limit, Eq. (l22l) implies that lepton number is locally conserved: L + R = 0. The Higgs density H, created 
by the CP-violating source, is converted into L through tau Yukawa interactions, until chemical equilibrium 
(I2TI ) is reached, when 

Case (ii): take Dr 00, keeping Dl finite. Any R density created by tau Yukawa interactions instantly 
diffuses away to 2; = ±00; therefore, we set R = 0. Now, tau Yukawa chemical equilibrium ((2TI) implies 

Liz) = ^H{z). (26) 
kn 

In other words, tau Yukawa interactions will enforce chemical equilibrium locally. Since the RH lepton 
density is diffusing away, reducing the local R, more conversion of H into R and L occurs to compensate, 
thereby resulting in more LH lepton density. This conversion ceases when Eq. ( |26l ) is reached. Therefore, a 
large diffusion constant for RH leptons enhances the density for LH leptons. This enhancement, maximized 
for Dr 00, is at most a factor of 

^^-^ ~ 0(5) . (27) 
kR 
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Both cases agree with Eqs. (l24l ). setting AL, Ai? —>■ 0. 

Next, consider the case of physical relevance, where D^j ^ Dl, but keeping both Dr, Dl finite. Close 
to the bubble wall, LH lepton density will be enhanced, as argued above. However, far from the bubble wall, 
an additional effect occurs: RH lepton density, having diffused far into the unbroken phase, is converted into 
L and H by tau Yukawa interactions. This effect suppresses L. Close to the bubble wall, Higgsinos created 
by the CP- violating source (H > 0) will be converted into LH leptons (L > 0) and RH anti-leptons {R < 0), 
and then, far from the wall, the RH anti-leptons will be converted into LH anti-leptons, thereby suppressing 
L. This physics is incorporated in the non-local corrections AL and Ai?, which we now consider. Using 
Eqs. ((2n[22l[24)) . we can derive differential equations for these densities: 

- Dlr AL" + AL' = , ^2 ^ V H' (28a) 

KHikL + KrY DlR 

- Dlr AR" + AR' = , Dl-Dr ^, ^ ^^Sb) 

kH[kL + kRY Dlr 

where Dlr = {Dl Ul + Dr kR)/{kL + kR). With the boundary conditions AL(ibcxD) = AR{-^oo) = 0, 
the solutions to these equations are 

A«.) - ^ r MW) e-'-^'""- . (29b, 

These terms describe how regions of net lepton number can develop when Dr ^ Dl- Using Eqs. (I24I29I) . 
it is straight-forward to show that these solutions for L and R satisfy 



/oo 
dz{L + R) = {) 
-oo 



(30) 



conserving lepton number. 

In our numerical study, we find that the impact from AL and Ai? on the analytic computation of ur/s 
(given below) is only 0(10%). Since there are much larger uncertainties in the analytic computation, it is 
safe to neglect the non-local terms AL and Ai? from Eqs. (I24l) . 



B. Quark charge densities 



When top and bottom Yukawa interactions are in chemical equilibrium, the relations 
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are satisfied; cf. Eqns. dSllHl. Tliese equations imply that 

Q T B 

kq kt Kb 

First and second generation quark densities only couple to third generation densities, via strong sphaleron 
interactions, through the linear combination {2Q/kQ — T/kr — B/ks), as can be seen from Eqns. (fT4l) . 
Since this combination vanishes, third generation quark densities do not induce lst/2nd generation quark 
densities. Mathematically, if we impose Eq. 02] ). the Qi Boltzmann equation (|14c| ) becomes 

Vn^Q[-DgQ'l a -r,,Qi , (33) 

which, with the boundary conditions Qi{±oo) = 0, implies Qi{z) = 0. According to Eq. (fTTI) . we have 
Ui = Di = —Qi/2 = 0, for i = 1,2. Therefore, we conclude that all first and second generation quark and 
squark charge densities vanish in the presence of fast top and bottom Yukawa interations. Strong sphalerons 
only induce first and second generation densities in order to wash out an asymmetry between left- and right- 
handed quark chemical potentials; when bottom Yukawas are active, this asymmetry vanishes and strong 
sphalerons have no effect. 
Eqns. (I13I31I ) imply 

rr - rr 2kB + kQ 

T = ktH = - — ; — H 

kn ks + kQ + kr 

Q = kqH = — — — — — H (34) 

kn kB + kQ + kr 

kn kB + kg + kr 
The contribution to from third generation LH quarks is 

kg kB-kr , , 

n,,3 + n j3 = — ^ ; — H , (35) 

'^^ kH kB + kg + kr 

while that from first and second generation LH quarks vanishes. Let us contrast these results to previous 



work that neglected bottom Yukawa interactions 11411 



kq kB - 9kT 

+ „^ ^ ^ '->'^^'"' - + '"-^f- + ^^«) H. i = 1. 2 . (36b) 

"■f- "'^ kn kB + 9kQ + 9kT 

The formulae are completely different. Whereas in previous work significant baryon asymmetry could 
arise from first and second generation LH quarks, the presence of bottom Yukawa interactions completely 
changes the picture: no first and second generation quark density is created. In addition, with fast bottom 
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Yukawa interactions, the third generation quark charge vanishes when kx — ks, or equivalently m^^ ~ 

mr ; without them, this cancellation never occurs. 

Or 

Let us explain the physical origin of this cancellation. Suppose that the CP-violating source creates pos- 
itive Higgs/Higgsino density, such that H > 0. Due to hypercharge conservation, top Yukawa interactions 
convert Higgsinos and Higgs bosons into LH quark and squark antiparticles (driving Q < 0) and RH top 
quark and squark particles (driving T > 0), while bottom Yukawa interactions drive Q > and B < 0. 
Which effect wins is determined by whether T or B has more degrees of freedom available, according to the 
equipartition theorem. This is determined by the statistical weights kx and ks, which are governed by the 
masses mj^ and "i^^. When the masses are equal, we have ks — kx, suppressing n^s +n^3^. Similarly, the 
sign of (n„3 + n^3 )/H is positive or negative, depending on whether m^^ /"^b^ greater or less than unity, 
respectively. In the lepton-mediated scenario, we suppress the quark contribution by choosing mj^ ~ ^bji' 
In scenarios beyond the MSSM, it is also suppressed for mj^, m-^^ ^ T. 

C. Solving the Boltzmann equation 



In terms of H, the left-handed fermion charge density (1181) becomes 

nL{z) = -r- 1 ; r- ^i^) + i — i — f; 1 — ^i^) + i— ^L[z) , (37) 

kn kB + kq + kr ku kiDL + kRDR kL 

where AL is given in Eq. (I24al ). The first term is the contribution to from third generation quarks, 
while the second and third terms are contributions from third generation leptons. The lepton contribution is 
predominantly given by the second term only; the third term, as discussed above, is suppressed for v^, <C 1. 
This equation is the main result of this paper; from it, we infer several conclusions: 

• The lepton contribution is enhanced for m^^ <^ "m^f^, when kR is largest and ki smallest; {cf. 
Eqs. (I4I17I )). It is also enhanced for Dr » Dl. Its sign is fixed with respect to H, which in turn 
is fixed by the sign of the CP-violating source, as we show below. Therefore, in a lepton-mediated 
EWB scenario, where ni is predominantly leptonic, the sign of the CP-violating phase most relevant 
for EWB uniquely fixes the sign of n^/s, in contrast with the quark-mediated scenarios. 

• Left-handed charge arises from third generation quarks and leptons, and not first and second gener- 
ation quarks and leptons. The form of iil is qualitatively different than in previous treatments that 
neglected Tyi) and Tyr, where left-handed charge came from quarks of all generations, and not from 
leptons. 
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Furthermore, the quark contribution to vanishes for m^^ = "^^^> since ks = kr- Its sign is 
opposite to that of the leptonic contribution for m-^^ < and the same for mj^ > . 



We explore these implications in more detail numerically in SecHVl 

We emphasize that our conclusions are quite general, although it appears that our Boltzmann equations 
(11411161 ) have been specialized to the MSSM. In any extention of the MSSM, Eq. (|37] ) and its conclusions 
remain valid if the following conditions hold: (i) third generation Yukawa interaction rates are faster than 
the diffusion rate, and (ii) CP-violation is communicated to the first and second generation quark sectors 
solely through strong sphalerons. 

Since ns/s is determined by n^, all that remains is to solve for the Higgs charge density H. We 
can reduce the Boltzmann equations (11411161 ) into a single equation for H by taking the appropriate linear 
combination of equations 

([T4al) + 2 X (fT4b]) + ([B]) + , (38) 

such that the Yukawa and strong sphaleron rates all cancel, and expressing the densities L,Q,T in terms 
of H using Eqs. (I24I34I) . This master Boltzmann equation equation is an integro-differential equation for 
H{z), due to the presence of the AL term. Therefore, for simplicity, we treat AL perturbatively: first, 
we neglect AL in our solution for H, and then, given our solution H, we include the AL contribution in 
Eq. (|37] ) for ul. Neglecting AL, the master Boltzmann equation is 

Vn,H'-DH"= - FH + S, (39) 

where 

D = Dh + DQjKT - Kb) + Dl kl ^^^^^ 

1 + KT - KB + KL 

Y ~\~ Tmi "i" Trnfe "I" ^mr (40b) 

knil + KT - KB + Kl) 

5'<3'v _j_ ^cpv _ ^cpv _ ^cpv 
S = ^ — . (40c) 

I + KT - KB + KL 



1611 . the form of 



Although the expressions in Eq. ( [39l ) are identical to that in the established literature il4l 
Eqs. (l40l) is dramatically different. We note that there is no dependence on the first/second generation quark 
sector, owing to the fact that they do not participate in the dynamics which determines ni- 

To solve Eq. (l39l) analytically, we follow Ref. making the approximations (a) that the true spatial 



dependence of the chiral relaxation rates may be replaced by a step-function, so that we may write r{z) 
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r 9{z); and (b) that S{z) ~ for z < —Lw/2. For the symmetric phase, where z < —L^jl, we obtain 



(41) 



where 



A 



p-7+y 
^7+ 



J~Lw/2 



dy S{y) 



7- 



+ 



(42) 



Furthermore, we have defined 



7± = 



1 
2D 



± Jvl + ATD 



(43) 



We reiterate that although the form of Eqns. (I41H43I ) is similar to that in previous work 11411 . our results for 
D, r, and S are different, due to the modified structure of the Boltzmann equations in the presence of fast 
bottom and tau Yukawa rates. 

We now ask: was it safe to neglect AL in solving for HI Substituting our solution for H into Eq. ([29]) . 
we find that /S.L/H oc 1/Dr, in the limit that Dr ^ cxd. In short, in the physical limit where large RH 
lepton diffusion has the biggest impact upon ni, our solution is most accurate. There may be situations 
in general in which the impact of non-local corrections is not suppressed; we show how the Boltzmann 
equations may be solved in this case in future work 14411 . 



IV. LEPTON-MEDIATED ELECTROWEAK BARYOGENESIS: NUMERICAL RESULTS 

We now consider an MSSM scenario that illustrates some of the novel features discussed in Sec. [Till 
As we will see, the picture here is that the BAU is induced predominantly by leptonic left-handed charge: 
hence, lepton-mediated. The key parameters that govern the behavior of this scenario are (i) tan /3 > 20 and 
pseudoscalar Higgs mass (at zero temperature) < 500 GeV, ensuring Tyr, Tyr <^ Tdiff, and (ii) right- 
handed top and bottom squarks with approximately equal mass, thereby suppressing the quark contribution 
to riL- Here, we take both squarks to be light, with 0(100 GeV) masses, since a strong first order phase 
transition requires a light top squark. 

Although we work within the context of the MSSM, many of our conclusions are much more general. 
In EWB scenarios beyond the MSSM, light squarks are not required for a strong first order phase transition 
(see e.g. Refs. i4d. l4lL 14211 ). Even if the squarks are very heavy, EWB is still mediated by leptons occurs as 
long as the previous two conditions are met. 

In this section, we first summarize the parameters of the lepton-mediated EWB scenario. Next, we 
compute the bottom and tau Yukawa interaction rates Fy^ and Ty^, showing for what regions of parameter 



16 





120 GeV 


Ml 


-(60 GeV)2 


T 


100 GeV 




Dq 


6/r 


Ml 


120 GeV 


Ml 


(100 GeV)2 


v{T) 


125 GeV 






lOO/T 


M2 


250 GeV 


Ml 


(300 GeV)2 


A/3 


0.015 




Dr 


380/r 


tan j3 


20 




150 GeV 




0.05 




Lw 


25/r 



TABLE I: Important parameters for lepton-mediated EWB scenario. 



space they are fast compared to Tdiff. Last, we numerically solve the system of Boltzmann equations (|14H16l l 
and compute the left-handed fermion charge density that generates n^/s. Our main result is Fig. (3) it 
illustrates how ni arises from leptons instead of quarks, how our analytic and numerical results agree, and 
how this scenario differs dramatically from previous work neglecting Tyi^ and Ty-r. 



A. Input parameters 

The computation of ub/s relies upon many numerical inputs, some described here and others described 
in the Appendix. We have evaluated the masses of particles during the EWPT assuming that electroweak 
symmetry is unbroken. This approximation is motivated by the fact that most of the charge transport dynam- 
ics take place outside the bubble in the region of unbroken symmetry. These masses receive contributions 
from the mass parameters in Table U and from finite temperature corrections, listed in the Appendix. The 
right-handed stop, sbottom, and stau SUSY-breaking mass-squared parameters are M|,, M^, and M^, re- 
spectively. The RH stop is required to be light to achieve a strong first order phase transition IsS^ ; taking 
the RH sbottom and stau to be light as well ensures that the quark contribution to rij^ is suppressed, while 
the lepton contribution is enhanced, in accord with Eq. (l37l) . We take all other squark and slepton (mass)^ 
parameters to be 10 TeV. 

For Higgs bosons, the story is more complicated. Again, we study the degrees of freedom assuming 
unbroken electroweak symmetry. The mass term in the Lagrangian is 



C^-(Hi^,H^]{ II I , (44) 




ml + l/ip + 




and the same for {H^, H^J) but with 6 — > — 6. The finite temperature corrections that restore electroweak 
symmetry are given by (see Table HII) 

= {l9l + \9l + \y^)T' (45a) 

.2 , 1 „2 I 1 „,2 , 1 „.2 \ rr2 



= ( -9'^ + -gi + -yt + -y',]T' . (45b) 
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(In the high T limit, there are no off-diagonal thermal corrections, since these corrections are proportional 
to dimensionful parameters.) We can re-express this mass matrix using the minimization conditions for 



electroweak symmetry breaking at T = 1431] : 



+ = m^cos^/3 + ^m|cos2/3 ~ -^m| (46a) 

'^d + I^P = sin^ /3 — cos 2/3 ~ + - m| (46b) 
h = sin I3q cos /? ~ , (46c) 

where niz and are the Z and pseudoscalar Higgs boson masses at T = 0. The approximations in 
Eqs. ( |46l ) follow assuming tan /? S> 1. Therefore, in this limit, the Higgs boson mass matrix (l44l ) is 
diagonal, with eigenvalues 

"^k = -\^l + ^u (47a) 
"^L = rn\ + ^m% + 5d ■ (47b) 

These are the Higgs boson masses during the EWPT. We note that if niA ~ 0(100 GeV), then H^, is also 
Ught. 

The top, bottom, and tau Yukawa interactions are proportional to the corresponding Yukawa couplings. 
At tree-level, these couplings are determined by 

iT^T nib mt 
Vr = ^ , Vh = ^ , yt = — —5 , (48) 

V cos p V cos p V sm p 

where v ~ 174 GeV is the Higgs vev at T = 0. However, quantum corrections lead to two complications. 
First, we include the QCD (QED) running of (y^) from the scale where its mass is measured to the 
electroweak scale mz', this reduces by a factor r/;, ~ 1.4 and has negligible impact on y,- [31]. Second, 
we allow for the possibility that is smaller than expected at tree-level, due to receiving large 
one-loop corrections enhanced by tan/3, denoted as 5i, and 5^, for which we include only the dominant 
contributions jsi . 
scale mz are 

yrimz) = ^.-1 , . , ^ , yb{mz) = ^.-1 , . , ^ • (49) 

V cos p [1 + Or tan p) v cos p [1 + Ob tan p) 

For parameters given in TableHl we find yb{rnz) — 0.33 and yAjuz) — 0.20. 



3211 . Including both of these effects, the Yukawa couplings evaluated at the electroweak 



2311 : the fact that Dr » Dl enhances 



The diffusion constants Di have been computed in Ref. 11221. 
the left-handed lepton charge, as discussed in Sec. Hill The bubble wall velocity v^, thickness L^, profile 
parameters A/3 and v{T) describe the dynamics of the expanding bubbles during the EWPT, at temperature 



18 



T i33h . The spacetime-dependent vevs are approximated by 

1 



v{z) 
f3{z) 



v{T) 



2 



1 — tanh 



3z 

Lin 



1 



A/3 



3]. 



1 + tanh 



3z 



(50) 
(51) 



following Ref. lit 

We consider a CP-violating source 59"^ arising solely from Higgsino-Bino mixing, enhanced for ^ = 

nn 

Ml, and calculated following Refs. 11161. 12511 : the relevant CP-violating phase is 4>i = arg(Mi/x), and the 
Higgsino and Bino thermal widths are ~ 0.025 T and ~ 2 x 10^^ T (assuming the Bino width is 
dominated by quark-squark loops) 



c-cpv ^ 



Numerically, we find 



0.24 GeV X (3'{z)v{zf sin 01 . 



(52) 



This "Bino-driven" CP- violating source requires |^| > 1/20, which is compatible with EDM searches 
provided we allow for arg(M2/u) / arg(Mi/i) [37]. The magnitude of S^~^ — and thus ub/s — is 
proportional to A/3, which itself goes as A/3 oc l/m\. Therefore, within this computation, viable EWB 
requires to be sufficiently light^ ; in the next section, we will see that light niA also increases the size of 
Tyh and Tyr. 

The CP-conserving relaxation rates wash-out CP-violating asymmetries within the broken phase. Com- 



puted following Ref. jlol . these rates are 



Th{z) ~ 3.8 X 10"^ GeV"^ x v{zf 
Tmt{z) ~ 3.0 X 10"^ GeV"^ x v{zf sin^ [3{z) 



r^,(z)~ (^^J cot' l3{z)T^t{z) . 
We neglect additional CP-violating relaxation rates from squarks, (s)leptons, and Higgs scalars. 



(53) 
(54) 

(55) 



B. Determination of bottom/tau Yukawa rates 



The leading contributions to Tyi, and Tyr arise from absorption and decay processes in the thermal 
plasma, shown in Fig.[T^. In the lepton-mediated scenario, the dominant processes are 



(56) 



' The computation of the CP-violating source is the subject of ongoing scrutiny. In other computations, there exist contributions 
to 5'^''^ that are not suppressed as niA ~^ oo (s^l . 
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FIG. 2: Contour plot of T^is/Tyb (left) and Tdiff/rj^r (right) in tan parameter space. Values of Tyb and Tyr 
include contributions from Hd qLbR and Hd ^ ^ltr only. Tau Yukawa chemical equilibrium is maintained for 
Tdiff / Tyr ^10; similarly for the bottom Yukawa. 

where ql = tL, ^L, and £l = tl, Vr- (Processes with left-handed squarks are suppressed since Mq = (10 
TeV)^; in addition, the decays < — > H are kinematically forbidden.) We compute these thermally- 



averaged decay rates following Ref. [17]; they are 

r 



— — ^Firuq^, mHa, nib J 
12 



J'2 



(57a) 
(57b) 



We refer the reader to Ref. 11711 for the general form of Zp, which is the form used in our numerical analysis. 
For the case of a scalar decaying into two fermions, •s-^ Vi it is approximately given by 



Ipimi, nis, m2) 



T3 



/ \5/2 



4 V27rT/ 



1 



TTT-l + m2 



(58) 



This simplied form for Tp is obtained by assuming > mi + m2 ^ \mi — m2\, and taking Maxwell- 
Boltzmann statistics for these particles; it is valid at the 0(25%) level for > 2{mi + 1712). 

The key parameters that govern Tyr and F^f, are tan /? (since the Yukawa couplings oc tan /? for 
tan /? » 1) and the masses of the particles in the decays in Eq. (l56l ). From Eq. (1581 ). we see that when 
the decaying scalars and tr are light compared to T, the rates Tyj- and Tyi, are largest. In particular, 
as shown in Eq. (|47] ). the mass of Hd during the EWPT is related to the masses of the Z and pseudoscalar 
Higgs at T = 0. Not only does light niA enhance the CP- violating source, it enhances F^;, and F^,-. We 

yb 



illustrate this fact in Fig. |2l here, we show the time scales for these rates, Tyb = T^^ (left) and Tyr = Tyr 
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FIG. 3: Left-handed charge densities for leptons (left panel) and quarks (right panel) that generate ns/ s, for lepton- 
mediated scenario. Solid (dashed) curves are numerical (analytic) results, as function of distance z from bubble 
wall. Shaded region denotes broken electroweak symmetry. Dotted curves are numerical results obtained neglecting 
tau/bottom Yukawa interactions. The effect of these interactions is to suppress LH quark charge, while enhancing LH 
lepton charge, thereby flipping the sign of til and ub / s compared to previous computations. 

(right), compared to the diffusion time scale Tdiff = D /vlj, as a function of tan (3 and m^. We only include 
contributions to these rates from decay. Bottom Yukawa interactions are in chemical equilibrium for 
'Tdin/Tyh ^ 10, whcn tan/? > 5 and niA ^ 600 GeV; similarly, tau Yukawa interactions are in chemical 
equilibrium for Tdiff /ry,- ^10, when tan /3 > 20 and niA < 500 GeV. For the parameters chosen earlier, 
we have 

Tyb ~ 0.28 GeV , Tyr ~ 0.12 GeV , Vyt ~ 2.4 GeV , (59) 

with Vyt computed in Ref. [ 17]; for comparison, t^^^ ~ 0.005 GeV. 

Scattering contributions, shown in Fig. [T}5, also contribute to bottom and tau Yukawa rates. They are 
suppressed in comparison by as or a^, but become the dominant contribution when absorption and decay 
are kinematically forbidden. For completeness, we also note another class of interactions mediated by third 
generation Yukawa couplings: F-term-induced four-scalar interactions, shown in Fig.[TJ;. However, one can 
show that if all Vyi (for i = t,b,T) interactions are in chemical equilibrium, then chemical equilibrium is 
satisfied for these four-scalar interactions as well. In the present work, we consider only the contributions 
to Tyi, and Fy-j- that arise from absorption and decay processes. Therefore, our calculation for these rates is 
a lower bound; the extent to which they lead to chemical equilibrium can only be enhanced by the inclusion 
of scattering processes in Figs.[T}5,c. 
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C. Numerical Solution to Boltzmann Equations 

We now solve the Boltzmann equations (|14H16I ) numerically for the lepton-mediated EWB scenario. In 
Fig. [3l we show the left-handed fermion charge density n l that arises from leptons (left) and quarks (right), 
for maximal CP-violating phase = — 7r/2. Our numerical results are shown by the solid curves, plotted 
as a function of the distance z to the moving bubble wall. The region of broken electroweak symmetry 
(denoted u / 0) is for z > 0, while unbroken symmetry is for z < (denoted v = 0). As promised, ul 
is predominantly leptonic, while the quark contribution is suppressed. The dashed curves are our analytic 
results (plotted only for z < 0), obtained in Sec. |llll Our analytic and numerical results are in good 
agreement. However, close to the bubble wall, there is some disagreement between numerical and analytic 



lepton charge densities. For \z\ < y Z) Ty-r ~ 2 GeV, the lepton density has not had enough time to reach 
chemical equilibrium; here, our analytic approximation is breaking down, as discussed in Sec.|IIl 
The resulting baryon asymmetry is 

{8 X sin 01 {nB/s)^^^ Bottom/tau Yukawas included 

(60) 
— 14 X sin 01 (f^B/s)cMB Bottom/tau Yukawas neglected 

where {nB/s)cMB = 8.84 x 10^^^ is the central value obtained from the CMB Mi. That is, a CP-violating 
phase sin (pi ~ 1/8 is required to give the observed BAU in this scenario. However, if we had neglected 
bottom and tau Yukawa interactions, the picture completely flips, as shown by the dotted curves in Fig. [3] 
In this case, the lepton charge vanishes (left panel) and the quark charge becomes dominant (right panel), 
thereby flipping the sign of the required CP- violating phase: sin 01 ~ —1/14. In the lepton-mediated 
scenario, the impact of bottom and tau Yukawa interactions is dramatic. If electric dipole moment searches 
uncover new CP- violating phases, such as 0i, the inclusion Fy^ and Fyr will clearly be essential in testing 
the consistency of supersymmetric EWB scenarios. 



V. CONCLUSION 



In this work, we have studied how the transport dynamics involving the tau and bottom Yukawa inter- 
actions impact the generation of left-handed fermionic charge n l that biases baryon number production via 
electroweak sphalerons. Previous work has neglected these interactions. However, we showed that these 
interactions are not negligible and can have a dramatic impact upon iil and ub/ s. When tau and bottom 
Yukawa interactions are in chemical equilibrium, the following effects occur: 

• Significant third generation lepton charge is generated. In contrast, without Tyr, no left-handed 
lepton asymmetry is generated. We showed how this charge is enhanced by efficient diffusion of 
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right-handed leptons. 

• No first and second left-handed quark charge is generated, and strong sphaleron processes are unim- 
portant. In contrast, without Ty^, a significant fraction of left-handed charge comes from first and 
second generation quarks, and strong sphaleron processes become important. 

• Third generation left-handed quark charge is suppressed if the right-handed top and bottom squarks 
have equal mass, or if their masses are large compared to the temperature. In contrast, without Tyh, 
this suppression does not occur. 

In light of these differences, the inclusion of Tyb and Tyr can have a strong impact on the computation of 
the baryon asymmetry. 

To verify our analytic conclusions, we calculated the rates Fyb and Tyr, and solved the full system of 
Boltzmann equations numerically. We considered a "lepton-mediated electroweak baryogenesis scenario", 
where all three right-handed third generation scalars have 0(100 GeV) masses, and all other squarks and 
sleptons have 0{TeY) masses. Here, we showed that all three left-handed quark densities are suppressed, 
and the baryon asymmetry is induced primarily from left-handed lepton charge. The CP-violating phase 
required for the observed baryon asymmetry is positive. In contrast, if we had neglected Tyi, and TyT, 
we would have arrived at a completely different picture: the baryon asymmetry would have arisen from 
left-handed quark density, and the required CP-violating phase would have been negative. The baryon 
asymmetry in the this scenario is strongly impacted by the presence of Tyi, and Tyr. Although we have 
focused on a lepton-mediated scenario, the leptonic component of ni can be important even when the 
quark component is not suppressed, enhancing or suppressing the total BAU. We conclude, therefore, that 
bottom and tau Yukawa rates are relevant and must be included in supersymmetric electroweak baryogenesis 
computations, both in the MSSM and beyond. 
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APPENDIX A: NUMERICAL INPUTS 

Here we summarize other numerical inputs needed to compute ns/s. In Table Ull we list various finite 
temperature mass corrections that relevant for EWB. The numerical values of gauge couplings used are gi = 
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TABLE II: Thermal masses Sm?' for active particles in the plasma during EWPT. The different contributions arise 
from thermal loops involving: (A) SM fermions and gauge bosons only; (B) Higgsinos, Winos, Binos, and RH stops; 
and (C) RH sbottoms and RH staus. In the lepton-mediated scenario, all three contributions must summed. 



0.357,(72 = 0.652,(^3 = 1.23 lH]. We note that ff^ have different finite temperature contributions; here, 
we treat these degrees of freedom as a single Dirac fermion with plasma mass equal to the average of these 
contributions. 

There are additional parameters of minor relevance to the EWB computation. We choose third generation 
SUSY-breaking triscalar parameters to be (y^ A), for i = t,b, r, with a common A = 7 TeV. In addition, 
we choose gluino mass parameter M3 = 500 GeV. These parameters are relevant for computing 6b,T in 
Eq. (|48] ). In addition, with these parameters, we can compute the T = squark and slepton spectrum; the 
lightest, mostly RH, third generation scalars have masses 



mr ~ 
tl 



103 GeV , mr ~ 93 GeV , m=, ~ 297 GeV . 



(Al) 



All other squarks and sleptons have 0(10 TeV) masses. The computation of n^^ is effectively insensitive to 
these heavy masses, since these particles have decoupled from the plasma before the EWPT. 

Lastly, we describe how to obtain ub/s after solving the system of transport equations for nL{z) — 
either analytically or numerically, as described in the text. The BAU is given by 



2D 0X4 



dz nL{z) e' 



-A_z 



(A2) 
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with 

A± = ^ (v^ ± ^JvI + ^DqTZ^ , (A3) 
with 7^ ~ 2 X 10^^ GeV, the wash-out rate for inverse electroweak sphaleron transitions (see e.g. Ref. llal). 



Baryon number induced by nL{z) is washed out for z » |Al^| ~ 30 GeV^^; this is another reason to 
neglect the non-local contributions, which contribute to n/, far from the bubble wall (z = 0). Lastly, the 
entropy density is given by 

s = —g,sT\ (A4) 
45 

where 5*5 ~ 131, for the parameters given. 
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